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New non-linear, spatially periodic, long wavelength electrostatic modes of an electron fluid oscil- 
lating against a motionless ion fluid (Langmuir waves) are given, with viscous and resistive effects 
included. The cold plasma approximation is adopted, which requires the wavelength to be sufB- 
ciently large. The pertinent requirement valid for large amplitude waves is determined. The general 
non-linear solution of the continuity and momentum transfer equations for the electron fluid along 
with Poisson's equation is obtained in simple parametric form. It is shown that in all typical hy- 
drogen plasmas, the influence of plasma resistivity on the modes in question is negligible. Within 
the limitations of the solution found, the non-linear time evolution of any (periodic) initial electron 
number density proflle ne{x,t — 0) can be determined (examples). For the modes in question, an 
idealized model of a strictly cold and coUisionless plasma is shown to be applicable to any real 
plasma, provided that the wavelength A ^ Amin(^^o, Te), where no — const and Te are the equilib- 
rium values of the electron number density and electron temperature. Within this idealized model, 
the minimum of the initial electron density rig (a;min,t = 0) must be larger than half its equilibrium 
value, no/2. Otherwise, the corresponding maximum ne(a:niax, t = Tp/^), obtained after half a period 
of the plasma oscillation blows up. Relaxation of this restriction on n^{x,t — 0) as one decreases A, 
due to the increase of the electron viscosity effects, is examined in detail. Strong plasma viscosity 
is shown to change considerably the density proflle during the time evolution, e.g., by splitting the 
largest maximum in two. 

PACS numbers: 52.30.-q,51.20.-|-d 



I. INTRODUCTION 



In this paper, which is an extension and generalisation of an earlier Letter [1], we present new features of an 
important mode in a plasma, i.e., electron plasma waves, also called Langmuir oscillations or space charge waves. 
The non-linear (fluid) features of these waves will be examined in full detail, along with approximate inclusion of the 
dissipative effects due to electron viscosity and plasma resistivity. 

Electron plasma waves were first examined, both experimentally and theoretically, in 1929 by Tonks and Langmuir. 
A simple theory was given [: ], treating the plasma electrons as non-interacting particles oscillating against the mo- 
tionless ions. This simple picture, neglecting random thermal motions, led to the equation of motion for each electron 
to be that for a harmonic oscillator with angular frequency 

= cjpg = 47rnoe^/me, (1) 

where e and me are the electron charge and mass, and hq — const is the equilibrium value of the electron number 
density. A similar equation was obtained for the electric field E under the assumption that V x E = 0, i.e., that 
there is no magnetic field associated with the wave. This means that the electric current due to electron oscillations is 
fully compensated by the displacement current E/(47r). This fact and formula (1) defining the frequency of electron 
oscillations are the most essential features of Langmuir oscillations in the "cold plasma approximation" , neglecting 
thermal motions. 

The infiuence of random thermal motions on small oscillations of electrons can be examined within the kinetic 
theory based on the linearised Vlasov equation. This was done correctly for the first time in 1946 by Landau [ ] who 
arrived at the dispersion relation in the form 

^2 ^ c^2^ + k^{ST,/m,) = w^Ji + 3{kXD)% (2) 

valid in the long wavelength limit, (kXjj)^ <C 1, where k = 27r/A is the wavenumber, A is the wavelength, A^ = Vth/^pe 
is the Debye length, vth = \/Te/me is the electron thermal speed, and is the equilibrium electron temperature 
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in energy units. In this long wavelength limit, the phase velocity of the wave, Uph = oj/k « ojpe/k is much larger 
than vth_. Landau's analysis also indicated that the electron plasma waves in this limit undergo a small coUisionless 
damping with decrement given by 
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(3) 



These results, along with the contour integrals used in the derivations, turned out to be very important for further 
progress in the kinetic theory of various waves in plasmas, see e.g., [4]. The decrement given by (3) is nowadays 
referred to as Landau damping. The dispersion relation (2) is commonly associated with the names of Bohm and 
Gross who derived it independently three years later [l]. However, their derivation, not based on either Vlasov theory 
or any other systematic method, had no chance of generalisations or finding damping. 

The physics behind Landau damping is nowadays well understood. It appears for any particle distribution function 
which decreases as one increases the particle velocity in some vicinity of the phase velocity of the wave. In that case 
there are more particles that are a bit slower than the wave and are accelerated on account of the wave energy than 
those decelerated. Large Landau damping is obtained if the phase velocity of the wave is close to the particle thermal 
speed where the slope of the distribution function is large. 

The simplest model that can be used to examine waves in a plasma is the fluid model in which electrons and 
each kind of ions are treated as fluids. In the small amplitude limit which allows for linearisation of the macroscopic 
equations for — uq, Vg, etc., one is dealing with linear and homogeneous partial differential equations for these 
parameters and the electric and magnetic fields, see e.g., [(>]. Very fine effects, like Landau damping, are beyond the 
scope of the fluid description. However, it is much easier to include non-linear and dissipative effects, the presence of 
external fields in equilibrium or equilibrium inhomogeneities. 

If the equilibrium parameters are space independent, the partial differential equations in question have constant 
coefficients and are satisfied by plane waves, where n,, — uq, Vg, E etc., are proportional to exp[z(k • r — wi)]. In the 
case of Langmuir oscillations, the waves are longitudinal (E and Vg have the direction of the wave vector k, due to 
V X E = 0, i.e., k X E = 0)). One can always choose the x axis of the coordinate system to have the direction of 
k. With this choice, k ■ r = kx and the vectors Vg and E have x components only. The real part of such a solution 
is one-dimensional, periodic in x with wavelength A = 27r/fc. Its profile is purely cosinusoidal. More complicated 
periodic profiles with the same A can be obtained by adding higher harmonics, with k replaced by nk, n = 1,2,... 
In the linearised theory, they will also be solutions, oscillating with lo ~ ujpe- The aim of this paper is to describe the 
behaviour of these ID solutions as one increases the amplitudes so that the non-linear terms cannot be neglected, and 
also with dissipative effects included. 

For a purely cosinusoidal initial profile of — uq and no dissipation, this problem was solved formally by Dawson 
in 1959 by introducing a Lagrangian coordinate and in more detail by Davidson and Schram in 1968 [ ]. The solution 
had a somewhat unexpected initial amplitude limitation. This will follow from our results as a very special case. 

In the linearised theory one can superpose plane plasma waves propagating in various directions. Such higher 
dimensional solutions are beyond the scope of our fully non-linear ID analysis. They are not so important as the 
ID Langmuir waves that are importantly recently revisited in connection with laser-driven plasma-based electron 
accelerators. Such accelerators are capable of supporting fields even in excess of 100 GV/m [b]. 

While in the kinetic theory of Langmuir waves as a rule one assumes that ions are motionless, in the fluid model 
one can easily include motions of all plasma species, to check this assumption quantitatively. To close the set of fluid 
equations involving the continuity and momentum transfer equations for all species, one has to postulate a pressure- 
density relation. If the polytropic one is assumed (pa/n^ = const), the relevant equations for the longitudinal waves 
propagating along the x axis are as follows [ ] (a = e for electrons and 1, 2, . . . , A'^; for ions). 

Dispersion relation: 

1-E , ^T^"/' ,2 =0^ .1^^"^^. (4) 

Complex amplitudes of the macroscopic velocity Va versus electric fleld E: 



*cjmQ 1 - 7a(?;tha/wph)^ ' 

where qe = — e, qi = ZiC, Utha = \/Ta/rna, fph — w/A:, Uao and are the equilibrium values of the number density 
and temperature in energy units for the a species (ueo = no). The thermal corrections (subtracted from 1 in the 
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denominators) must be small as compared to unity, to avoid strong Landau damping. Neglecting them we obtain 

r Ni 



2 2 



P 



(6) 



V.i = -Zi Ve- (7) 

rrii 

Eqs. (6) and (7) justify the assumption of motionless ions {rrii — oo). Under this assumption, Eq. (4) simplifies to 

LJ^ =ujl^ + k^ij,Tjm,). (8) 

This dispersion relation following from the fluid model will coincide with that from the kinetic theory, Eq. (2), if we 
choose 7e = 3. This value is thus appropriate for the fluid description of Langmuir waves, which can be interpreted 
as a one-dimensional adiabatic compression of the electron gas, associated with these waves. 

II. BASIC EQUATIONS 

For X dependent electron plasma waves and a motionless ion fluid, the model equations describing the electron fluid 
are: the continuity and momentum transfer equations along with the adiabatic pressure-density relation with 7e = 3, 
and Poisson's equation, see e.g., [')] (Gaussian units), 

dUf, d , 



merie 



= 0, (9) 

(dVe , 5We\ 9/4 dv^\ ^ 



Ve PO 



(11) 

47re(no - Ue), (12) 



dE 
dx 



where the equilibrium electron temperature Tg is in energy units, Vf, is the electron viscosity coefficient, 77 is the plasma 
resistivity, and an ideal gas equation of state for the electron gas in equilibrium is assumed, po — noTe- 

We assume that the first term on the right hand side in Eq. (10) is negligible as compared to the second one, i.e., 
that the wave is driven by the electric field rather than by electron pressure (cold plasma approximation). For small 
amplitudes that allow for linearisation, that will be the case if the thermal correction in the dispersion relation (2) is 
negligible, i.e., 

3 {kXof = f ^) ^ < 1: where (13a) 



, , r, kvth kvth vth TpVth ^ 

kAo = 27r — = sa = = , (13b) 

A ojpg u A 

and Tp = 27r/a; is the period of plasma oscillations. This condition is fulfilled in any plasma, provided that the 
wavelength A is sufficiently large as compared to the Debye length. In that case, A is also large as compared to the 
distance travelled by an electron moving with its thermal speed during one period of plasma oscillation, and phase 
velocity of the wave is large as compared to the electron thermal speed. The non- linear counterpart of condition (13) 
will be derived in Sec. VI. 

Analytical formulas for the electron viscosity coefficient Ve and resistivity ry involve several approximations, see 
Sec. VI. Using the latest values of the physical constants the relevant formulas given in the first reference of [9] 
take the form (two component plasma of electrons and one kind of Z = 1 ions, i.e., protons, deuterons, or tritons): 

o(TJeVl)^/2 , , 

^^ = 0.73 ^ " =4.01318x10-^ ^ /' , 14a 
AV2^Xce^ Ac/10 ^ ' 
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where Ac is a slowly varying Coulomb logarithm which will be treated as a constant: 



,'23.4 - ln(ny' if T, < 50 eV, ^, 

= , , i/2^_ix ..^ . (15) 



'0 

25.356 - ln(no'" T-i) if > 50 eV. 



(Typically w 10-20). 

Neglecting the electron pressure term, we can write Eq. (10) in the form 



where ly = (4/3) i 



dVe ^ dVe e , ^ 9^^e e^rLQ 

at ax nif. mf.n^ ox'^ rrie 



III. NEW VARIABLES 



The nonlinearity on the left hand side of Eq. (16) can be eliminated by introducing Lagrangian coordinates: xq(x^ t), 
the initial position (at i = 0) of an electron fluid element which at time t was at x, and time in the electron fluid 
rest frame, r(= t). The basic transformation equations between Eulerian and Lagrangian coordinates are (see [ 1 1] for 
more detail) 

d d d r 

v„ — , X = xq + Ve{xo,T') dr' . (17) 



dr dt dx jQ 
The continuity equation (9) in the electron fluid rest frame is now simply 

dx 

rieTj — = Ueoixo) = ne{xo,t = 0). (18) 
OXo 

This indicates that the nonlinear operator n^^d/dx (= n'^Q d / Oxq) in Eq. (16) gets a bit simpler in the Lagrangian 
coordinates. However, to make it linear, a new space variable s has to be introduced satisfying: 



Id Id Id 



ne{x,t)dx neo{xo) dxQ nods 
The second equality here implies the definition of the auxiliary variable s: 



(19) 



s = n^^ neo{xo) dxo, (20) 



which turns out to be proportional to the integral of the initial electron density profile. 
Using Eqs. (17) and (19), Eq. (16) takes the form 



dvp e V d 

E 



'ne\ dve 
■.uq) ds 



V -Ve = 0. 21 



dr rrie ^Tief^o ds 

An important point is that E can be linearly expressed in terms of Ve- Indeed, using Eqs. (12) and (9) it follows that 

dE 

— = ineneVe, (22) 

which expresses the fact that the electric current, ^erieVe, is compensated by the displacement current {4TT)^^dE/dt. 
Adding this equation to Ve times Eq. (12) we end up with 

dE 

— = AnenoVe, (23) 
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where the right hand side is Unear in Ve as promised. Therefore, if we hnearise the viscous term in Eq. (21) by 
replacing Ue by no, differentiate this equation d/dr and use (23), we obtain 

OT^ nieno OS'' or in or ^ 

where is defined by Eq. (1). Linearisation of the viscous term is a simpUfication justified by the fact that this term 
only corrects the motion, whereas the driving electric force and the convective non-linearities in the fluid equations 
are taken into account exactly. Furthermore, the formula for the viscosity coefficient, Eq. (14a), is approximate, up 
to a factor of two or three. 

Eq. (24) is a linear partial differential equation (PDE) for We(s,T) with constant coefficients. Solutions of such 
PDEs are any superpositions of the normal modes exp[i(fcs — fir)], for which Eq. (24) leads to the dispersion relation 



n = ±Lo{k) - ia{k), 
a{k) 



2meJio 2me ' 



w(fc) = ^ujl,-a^{k). 



(25a) 
(25b) 

(25c) 



Assuming that k is real and superposing the normal modes corresponding to the plus and minus signs in VL given by 
Eq. (25) we obtain four real solutions: 



(26) 



where f,g = sin or cos. 

For each r, Vg given by Eq. (26) is a periodic function of s with wavelength A — 271 /k. By adding higher harmonics, 
obtained from Eq. (26) on replacing k — > nfc, n = 1, 2, . . ., and multiplying by an amplitude, any solution periodic in 
s with wavelength A can be obtained. Replacing asin(wT) + &cos(a;T) = Acos(a;r + 1^9), the solution in question can 
be written as 



g-a(nfc)r 1^^^ cj(nfc)T + iy9i„ sin(nfcs) + y42n COS uj{nk)T + if 2n cos(nfcs)|. 



where Ain, ipin, ^2n, and ip2n are arbitrary constants. 



IV. GENERAL SOLUTION IN PARAMETRIC FORM 



(27) 



Our equations and final results take a simple and universal form if we introduce dimensionless quantities: 

X — kx J s — ks^ t — ujipfDij^ T — ujipfDT ^ 
uj(nk) I ^ _ 2- - 

\ 1/2 

-10 



= 1.31642. ru - , , 

10 (re[eV])3/2- 
2j/efc2 fiTefeVl 



2meUJpe 



Aln 

no' 



A - 

A2n — ■ 

Wph 

Eek 



E : 



(28a) 
(28b) 

(28c) 

(28d) 

(28e) 

(28f) 
(28g) 
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Thus, Eqs. (12), (19) and (23) are now: 



dE 
dx 



= 1 



ds_ 
dx 



dE 
'df 



(29) 



Integrating the first two over dx and the last one over df and using (27), we end up with equations which define all 
relevant quantities: x, E, he and Ve in terms of s and f {=t). Dropping bars for simplicity, the final results for the 
dimensionless quantities become: 



x{s,t) 
E{s,t) 



s + E, 

OO 



Ain9in(t) sin(ns) + ^2«52n(0 cos(ns) 
Ainguiit) cos{ns) ~ A2ng2n{t) sin(ns) 



n=l 



C30 

We(s,t) = ^ e"""* ^Ainfin{t)sin{ns) + A2nf2nit)cos{ns) 



COs{uJnt + (fjn), 3 = 1,2, 

anfjnit) - ujn sin(cj„i + (fijn) = cos(a;„< + (pjn + Vo«), 
arctan(w„/Q!„). 



(30a) 
(30b) 

(30c) 

(30d) 



(30e) 
(30f) 



The identity in Eq. (30c) is due to the fact that af^ + uj^ = 1, and therefore one can always find such (pon that 
cos((^o«) = an, and sin(v?on) = ^^n- 

The independent parameters Ai„, ipin, ^2rn ^nd ip2n must be chosen so as to ensure reality of all dependent 
parameters defined by Eqs. (30), x, E, Ue, and v^. 

For modes with n not too large, such that a„ = n^i/ + 7] < I, ujn and ipQn are real and positive, and fjn{t) and 
gjn{t) are purely oscillating cosine functions, with period 2tt/uj„ and arguments shifted by fon- 

For higher modes, either a;„ = ipon = 0, if a„ = 1, leading to fjn{i) = djnit) = cos(ipjn) = const, or a;„ and (pon 
become pure imaginary (= and i|(^on|): if ctn > 1- Only the latter case requires further attention, whereas in 

the remaining cases (q;„ < 1) we can choose Ajn to be real, j = 1, 2. 

Thus if a„ > 1, real values of the phase shifts fjn in general lead to fjnit) and gjn{t) being complex and having 
time dependent phases. Therefore in that case one cannot produce real results by an appropriate choice of complex 
coefficients Ai„ and A2n- The only exceptions are ipjn = or — ipon and ipjn = (0 or — (pon) T 7r/2. 

If we choose 

= 0, the resulting fjnit) = cosh(|LLi„|t) and gjn{t) = cosh(|LLi„|t + \(pon\) will be real. A similar 
situation will arise if we choose (pjn = —y^on, for which fjn{t) = cosh(|a;„|i — |iy9o„|), gjn(t) = cosh(|u;ji|t). In all these 
cases we can choose Ajn to be real, j = 1, 2. 

If we choose (^j„ = Ti'/2 (and also for ipjn = —'fion T ""/S), it is convenient to replace the amplitudes Ajn by 
Ajn/ujn- With this choice Eq. (30) will hold if new definitions of the functions fjn{t) and gjn{t) are adopted, i.e., 



gjn{t) 



± 



± 



sin(cLi„t) 



sin(a;„i) 



- cos(a;„i) 



(31a) 
(31b) 



fjnit) 

gjn{t) 



± 



= ± 



sin(w„f) 



an 

UJ„ 

sin(w„t) 



cos(w„t). 



(32a) 
(32b) 



where must be replaced by t if a;„ = 0. It can be seen that the right hand sides of Eqs. (31) and (32) are well 

defined and real for any n (if a„ > 1, cos(w„t) = cosh(|w„|t) and sin(w„t)/a;„ = sinh(|w„|i)/|w„|). Their values at 
t = are either zero or ±1. Thus, if Eq. (31) holds for both j = 1 and j = 2, this corresponds to the unperturbed 
value of the initial electron velocity {ve{x,t = 0) = 0), and maximally perturbed initial values of the electric field E 
and rig. If in that case {(pjn = ^t^/^) we choose A2n = and Ain = An, we will reproduce our earlier results given 
in [1]. Similarly, Eq. (32) corresponds to the unperturbed initial values of the electric field and the electron density 
{E{x, i = 0) = and ne{x, t = Q) = 1), and maximally perturbed initial value of the electron velocity. 
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Real values of fjn (t) and gjn (t) for a„ > 1 will also be obtained if we choose ipjn to be pure imaginary. 
Eq. (30c) is only meaningful if the sum is smaller than unity, which imposes a limitation on the amplitudes Ain 
and A2n- 

If ipin = ^2m the mode in question is a product of a function of t and a function of s, i.e., represents a standing 
wave. 

Note that all modes, both with periodic fjn{t) and gjnit) (for a„ < 1) and aperiodic ones (for q;„ > 1) are damped 
exponentially as t — >■ oo (with decrement a„ if a„ < 1 or Q!„ — — 1 = [-^/a^j — 1 + an]^^ if a„ > 1). 



V. SOLUTION GIVEN IN TERMS OF PHYSICAL VARIABLES 



One can eliminate the parameter s from Eqs. (30) thereby making x and t the independent variables. This parameter 
has to be determined in terms of x and t from Eq. (30a) and used in the remaining Eqs. (30). While numerically 
this is a simple task, analytical formulas are complicated, see Eqs. (33)-(37). At the same time, the parametric form 
(30), which involves simple elementary functions, can also be used to plot E, n,, and Ve as functions of the physical 
variables x and t (see Figs. 1-5 and 8-11 obtained by using "ParametricPlot3D" of Mathematica 7). 

The electron density ne{x, t) is a periodic function of x with wavelength 27r. It can be expanded in a Fourier series 
in X with time dependent Fourier coefficients: 



where (jie dx — ds) 



and similarly, 



n^ix, t) = 1 + |^-Bim(0 cos(mx) + B2m{t) sin(ma;) 

m—l 



Blm{t) 



2tt 



ne{x, t) cos{'mx) dx 



2tc 



|m s + E{s,t) I 



ds, 



B2m{t) 



1 . r r 11 

— / sin-^ m s-\-E{s,t) J- ds. 



(33) 



(34) 



(35) 



with E[s,t) given by Eq. (30b). 

Equation (33) along with the first equation in (29) integrated over dx leads to the Fourier expansion of E{x,t): 



E{x,t) = V — \~Bi^{t)sm{mx) + B2m{t) 



cos(rna;) 



(36) 



And finally, Eqs. (22) and (28) lead to 



ne{x,t)Ve{x,t) = 



dE{x,t) 
dt ■ 



(37) 



This equation along with (3G) and (33) defines Ve{x^t). 

The freedom in the choice of the constants Ain, fin, A2n, and (p2n can be reduced from four to two if we prescribe 
the initial electron density profile neo(a;) = n^{x,t = 0). Putting t = Q in Eq. (30c) and using standard formulas for 
Fourier coefficients we arrive at {n~'^ ds = dx) 



^l«5ln(0) 



— / cos|n[a; — _E(x, 0)] I da; 

ITT Jq L J 

n I neo[x') dx' 



nil 
1 

nTT 



cos 



dx. 



(38) 



where E{x,Q) is given by Eq. (36) and we have chosen s(0,0) = 0. An analogous formula for A2„S'2n(0) can be 
obtained from Eq. (38) if we drop the minus signs in front of the integrals dx and replace the cosine functions by 
sines. 

In general the integrals in Eqs. (34), (35) and (38) (first line) are not expressible in terms of elementary functions, 
but if the sum over m or n is truncated at some M or N , one can calculate as many integrals as needed (numerically). 
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If either M = 1 or iV = 1 (first harmonics only), and additionally -821(0) = or A21 = (i.e., neo{x) or ^is,t) 
as a function of s are even functions) the relevant coefficients are expressible in terms of Bessel functions. Thus using 
the identity [ ] (p. 423) 

/ cos{ax — zsinx) dx — TrJa{z), (39) 
Jo 

we obtain 

AmginiO) = -(-l)"+V„(nBn(0)), (40) 
n 

and = 0, if i3i,„ (0) = for m > 1 and = 0. 
Similarly, 

Bi„(t) = 2J™[m^ne-"^*5iiW], (41) 
and B2m{t) = 0, if Ain = for n > 1 and A2n = 0. 



VI. APPLICABILITY CONDITIONS 



As for applicability of Eqs.(14a), (14b), one assumes that the plasma is quasineutral (rtg « rii) and the distribution 
functions are not far from local Maxwellians. Other approximations come from standard assumptions of the Chapman- 
Enskog method. Therefore, uncertainty factors of two or three cannot be excluded [: ]. 

If the requirement that the electron pressure term in Eq. (10) is negligible as compared to the electric force is 
expressed in terms of the dimensionless quantities (28), it takes the form 

'I^oT' — fc ' ^ l)x^'^''lh' 
in view of (29). Differentiating Eq. (30c) d/ds and comparing the result with Eq. (30b) it can be seen that 

where both sides differ in higher space harmonics only. Finally the requirement in question takes the form 

which is the non-linear counterpart of condition (13). In the linearised theory, where the dimensionless Ug ~ 1, this 
condition takes the form of Eq. (13). Expressing Tg in eV and calculating Te[eV]/(noA^) from Eq. (28d), we end up 
with 

«emaxl''7« 0.04. (43) 

This indicates that at least one of the coefficients v ox r] must be much smaller than 0.2 (jiemax > 1)- Note that for 
given plasma parameters uq and Tg, is defined by Eq. (28c) and v is related to A by Eq. (28d), i.e., either v oi \ 
can be prescribed. 

In Table I we present 77 and ^/uX as functions of ng and Tg. As oscillating results are obtained only for v < 1^ 
values of the wavelength A must be greater than Amin = X{v = 1) given in the last column of Table I. Values of A in 
cm can be obtained by multiplying the last column of Table I by V v~^. For example, for v = 10~^", V = 10^ and 
the last column of Table I gives the values of A in km. Similarly for v = 10^*, it gives A in m, etc. 

Table I indicates that for realistic plasmas, i] is always much smaller than one. Therefore the requirement (43) 
is nearly always fulfilled. Thus, if ly (< 1) is close to unity so that the viscous term in Eq. (21) is significant, both 
quasi-neutrality and linearisation of Eq. (21) require that the dimensionless ne{x,t) cannot be significantly different 
from unity («e„iax ~ !)• the opposite limit of ^ 1, the viscous term in Eq. (21) is insignificant and there is no 
need to worry about its linearisation or quasi-neutrality. In that case Ug max can be large, but condition Eq. (43) will 
be fulfilled for reasonable values of riemax (if "^emax^ < !)■ 
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TABLE I: Dimensionless plasma resistivity rj and square root of dimensionless viscosity times the wavelength y^A in some 
typical plasmas. Choosing some value of u (< 1) and multiplying the last column by \fv~^ one obtains A in cm. 



Plasma Type 


no cm ^ Te eV 


Xd cm 7] 


Amin = A cm 


Interstellar plasma 


1 1 


7 X 10^ 3.08 X 10~^° 


9.37 X 10^ 


Gaseous nebula 


10^ 1 


20 8.30 X 10"'' 


5.71 X 10^ 


Solar corona 


10*^ 10^ 


7 3.03 X 10"^° 


9.44 X 10^ 


Diffuse hot plasma 


10^^ 10^ 


7 X 10"^ 2.12 X 10"^ 


35.7 


Solar atmosphere, gas discharge 


10^^ 1 


7 X lO"'^ 9.59 X 10"* 


5.31 X 10"^ 


Warm plasma 


10^'' 10 


2 X 10"* 4.47 X 10"^ 


7.78 X 10"^ 


Hot plasma 


10'" 10^ 


7 X 10"* 1.82 X 10"^ 


1.22 


Thermonuclear plasma 


10^5 10* 


2 X 10"^ 7.20 X lO"'' 


61.3 


Theta pinch 


10^'' 10^ 


7 X lO"'^ 1.52 X 10"^ 


4.22 X 10"^ 


Dense hot plasma 


10^* 10^ 


7 X lO"'' 1.22 X 10"* 


1.49 X 10"^ 


Laser plasma 


10^° 10^ 


7 X lO"'^ 9.13 X 10~* 


5.44 X 10"^ 


If in a realistic plasma (where 
mated by an idealized model of a 
with Oi.fl = and io^i = 1 give an 


r] is always small) we also choose <C 1, the plasma behaviour will be well approxi- 
plasma being strictly cold and collisionless (Te = v = j] = 0). In that case, Eqs. (30) 
exact and general (periodic) solution of our non-linear initial equations (9)-(12): 


a;(s,f) = 


s + E, 




(44a) 


E{s,t) = 


oo 

[^iri cos(i + ipin) sin(?is) ^ 

71=1 


- A2n COs(t + (p2n) COs(ns)] , 


(44b) 




oo 

1 + ^ ?^[^lr^ COs(i + (fin) COS 
n=l 


;(ns) - cos(i + (p2n) sin(7is)] , 


(44c) 


i'e(s,i) = 


oo 

- ^ [Ain sm{t + ipin) sm{ns] 

n=l 


1 + A2n sin(t + ip2n) cos(ns)] . 


(44d) 


As E{s,t) in Eqs. (44a), (44b) is a periodic function of t, the Fourier coefficients Bim{t) and B2m{t) given by (34) 
and (35) are also periodic, with period tt. This in turn implies the same periodicity of all plasma parameters. They 
will all oscillate, with frequency (in real time) equal to the plasma frequency cjpg. The fact that this frequency is 
independent of the amplitude of nonlinear oscillations is a characteristic feature of an ideal plasma, i.e., of any real 



plasma for which v,rj ^ 1. This fact is not new but is demonstrated here for the general solution (44). 

Another characteristic feature of the general solution (44) is that the dimensionless ne{x,t) is always greater than ^ 
and blows up if n^ix, t) ^ somewhere. One obtains rie oo when the sum in Eq. (44c) tends to —1 (its unreachable 
minimum) while its value at time shifted by tt (the half period) tends to 1 (its unreachable maximum). The fact that 
in a cold and collisionless plasma nemin must be greater than ^ must be kept in mind when prescribing the initial 
electron density profile nf,o{x) = ne{x,t = 0). Any minimum in the initial density distribution tending to ^ produces 
an exploding maximum after a half period of plasma oscillation. However, when the plasma viscosity v is increased so 
as to introduce a noticeable damping, growth of the maximum in question is reduced, and the corresponding minimum 
in the initial distribution can fall below i. This will be illustrated in Section VII. 

Table I indicates that for a typical thermonuclear plasma with = 10 keV, i] is much smaller than one and so 
condition (43) for the plasma being "cold" is fulfilled if n-emax < 1- This would be so even if one could increase the 
number density ng of such a plasma to the value reached in a laser plasma, uq = 10^° cm"'^, which would result in 
77 = 1.52 X 10"^. 

VII. EXAMPLES 

In view of the fact that for realistic plasmas rj is always much smaller than one, the actual value of 77 will have no 
visible effect on the plots of ne{x,t), etc. where it will be taken equal to zero. The plots will depend on whether or 
not 1/ is much smaller than one. In any case, values of A in cm can be obtained by multiplying the last column of 
Table I by Vi^. 
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FIG. 1; (color online). Plot of ne{x,t) when only An — 0.65 is nonzero, and u — 0. This describes all situations with f <^ 1. 
The corresponding A, for given no and Te, can be found from Eqs. (28c), (28d) or Table I. 



In Figs. 1 and 2 we present typical examples of the electron density evolution for a small number of modes included 
in Eqs. (30). Various spatially periodic structures can be produced but the initial density distribution follows. 

We can prescribe the initial density distribution ne{x, t = 0) either directly, or by defining initial values of its Fourier 
coefficients Bim{0), -B2m(0). In the latter case periodicity is guaranteed and in any case, the Fourier coefficients Ai„ 
and can be calculated by using Eq. (38). In the simplest case of only the first harmonic present, Ain are expressible 
in terms of Bessel functions, see Eq. (40). Otherwise numerical integration is necessary which, however, is a simple 
numerical task. Examples are shown in Figs. 3-5. 

The particular solution (30) in which Ain is given by (40) with ipin = — 7r/2, reduces to that of [ ] if = 77 = 
0, though in a different notation. Denoting, as in [ ], the amplitude i?ii(0) of the initial density deviation from 
equilibrium by -Bi(O), the known condition Bi{0) < ^ follows from the identity Jn{nz) — 2(i-z) ['-] (P- 935). 




FIG. 2; (color online). Plot of ne{x,t) when only An — 0.4 and A12 — 0.2 are nonzero, and u — 0.02. u ^ — 7.07, see 
Table I.) The viscous damping is evident after a single period. 
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FIG. 3: (color online). Plot of n^{x, t) when only Bi(0) = Bii(O) = 0.55 is nonzero, v = 0.015 and = 40. = 8.16, see 

Table I.) Note that in the presence of even weak viscosity -Bi(O) can exceed 1/2, see Fig. 6. 



For -Bi(O) = ^, ne{x,t — 0) reaches its minimum equal ^ a,t x — tt and ^ becomes zero a,t s — t — n (ue — > oo). 
This illustrates the general property of an ideal plasma formulated at the end of Sec. VI. 




FIG. 4: (color online). Plot of neix,t) when only _Bi(0) = Bii(O) = 0.55 is nonzero, ly = OA and A'' = 40. {ViP^ = 3.16, see 
Table 1.) The observed bifurcation of the maximum under strong viscosity is a new nonlinear effect. 

The behavior of the solution in question for > is shown in Figs. 3 and 4. In a viscous plasma Bi{0) can 
exceed 1/2, see Fig. 6. This figure assumes <C 1 and gives the minimal admissible value of v for given Bi(0). For 
I' = ^'min [^i(O)] , the smallest minimum of n~^{x = TT,t) a.t t = tmin becomes zero. Note also that if i/ is sufficiently 
large, see Fig. 4, where v = 0.1, a new nonlinear effect can be noticed, i.e., the largest density maximum splits in two, 
with a saddle point between the peaks. This effect is due to the presence of the integer in front of in Eq. (28b). 

The integral dx in Eq. (38) is elementary if nf,o{x) is a sequence of step functions. An example of this is shown 
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FIG. 5: (color online). Plot of ne{x,t) when only Bii(O) = 0.4 and i3i2(0) = 0.3 are nonzero, u = 0.015 and A*' = 40. 
= 8.16, see Table I.) 




FIG. 6: Minimal admissible value of v as a function of -Bi(O). 



in Fig. 7, where rig^Q^x) — 1 is given in the form of two pulses, pi and P2, emerging from a reference level r. The 
parameters pi, p2, and r can be positive, negative or zero but none of them can be smaller than —1, as that would 

«eoCx)-l 



0.5 
0.4 
0.3 
0.2 
0.1 


Pi 

P2 

X 


1 X 


2 ^ 


"3 J 


■4 


0.1 
0.2 


r 




I 3 


4 ; 


i 6 



FIG. 7: An example of possible initial electron density deviation from equilibrium, ne{x) — 1. Here: xi = 1, X2 = 2, xs = 3.5, 
X4 — 5, pi = 0.5, P2 = 0.3, and r = —0.251 as calculated from (46). 
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FIG. 9: (color online). Plot of ne{x,t) for a two pulse profile, pi = 0.5, p2 = —0.4, u = and TV — 100. 
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make < 0. Using Eq. (38) (the second line) we easily find: 

Aingin{0) = ^-irj-^— [sin(7isi) + sin(ns3) - sin(ns2) - sin(ns4)] 

nn^ L 1 + r 

H |"sin(ns2) ~ sin(nsi)l H |"sin(riS4) — sin(ns3)l \ . (45) 

1+Pi 1+P2 J ' 

where 

si = {l + r)xi, 

52 = Si + {l+pi)ix2~Xi), 

53 = S2 + {1 + r)(.X3 - X2), 
Si = S3 + {1+P2){X4,-X3). 

^2nff2n(0) wiU be given by the right hand side of Eq. (45) if we replace sin by cos. 

The integral Jq^ {neo — 1) dx must be zero, leading to a constraint on the independent parameters. Solving it for r 
we obtain 

^_ Pl{x2 - Xi) + P2{X4: - X3) ^^g^ 
2tT - {x2 - Xi) - (X4 - 0:3) 

Examples are shown in Figs. 8 and 9. 
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If the integral dx' in Eq. (38) is expressible analytically, the integral dx in this equation can easily be found 
numerically. For example, for the saw-tooth initial density profile we obtain (ipin = — 7r/2): 

Ain = — —[ cos\nx\(l + b-bx/Tr)]\ dx, = 0, (47) 

nn Jq I ^ -'J 

where b is the "amplitude" of the initial density perturbation {rieaix = 0) = 1 + 6). The results are shown in Figs. 10 
and 11. 

Note that if we linearise Eq. (16), by neglecting Vedv^/dx on the left hand side and replacing n,, by no on the right 
hand side, the resulting equation will be identical with the linearised form of Eq. (21) if we replace s by a; and r by t. 
This means that V(.{x, t) for the linearised problem will be given by Eq. (27), or its dimensionless equivalent Eq. (30d), 
if we replace s by a; and r by t. Then, using the third equation in (29) (with r — > t) and the first equation we can 
find E{x,t) and ne(x,t). Finally, E{x,t) will be given by Eq. (30b) with s replaced by x, and for Ue we will get 

oo 

ne{x,t) = 1 + ^e"""*n[Ai„5i„(i)cos(na;) - A2ng2n{t) sm{nx)] . (48) 

n=l 

Putting here and in Eqs. (30b) and (30d) Ain — for n > 1, A2n = 0, and (pu = — 7r/2 we obtain 

E{x,t) = -yliiGi(i)sina;, (49a) 
ne{x,t) = 1 + AiiG'i(i)cosj:, (49b) 

Ve{x,t) = Ane-"i* """^sinx, (49c) 

UJl 

Gi{t) = e-"i*U^^^^^^+cos(wii)l, (49d) 
L wi J 

where wi = \/l — a\, ai = v + r]. This result becomes a linear counterpart of the nonlinear solution discussed earlier 
(only Bii{0) nonzero) if we put An — Bi{0). The amplitude An can take any value between —1 and 1, but there is 
no maximum amplification shown in Fig. 3. There is also no maximum splitting for large enough viscosity shown in 
Fig. 4. Both are clearly nonlinear effects, not present in the linearised theory. 




FIG. 11: (color online). Plot of ne{x,t) for a saw-tooth initial density profile, b — 0.55, ly — 0.1 and A*' — 40. {V ~ 3.16, 
see Table I.) Again the bifurcation of the maximum under strong viscosity can be seen. 
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VIII. CONCLUSIONS AND FINAL REMARKS 



The ID Langmuir waves are importantly recently revisited in connection with laser-driven plasma-based electron 
accelerators. A full non-linear fluid description of these waves, with the dissipative effects included, is given by 
Eqs. (30). They define the dynamics of the macroscopic parameters of the wave in a simple parametric form, only 
containing trigonometric and exponential functions. This simple form can be used to represent graphically the time 
evolution of the wave, ne{x,t), etc., while direct analytical formulas for these quantities, Eqs. (33)-(37), are much 
more complicated. 

Our analysis indicates that in real plasmas, both space and laboratory, the influence of plasma resistivity and 
electron pressure forces on the waves in question is negligible. 

The role of electron viscosity can be noticeable as one decreases the wavelength A of the Langmuir wave excited in 
a plasma. However, for A sufficiently large as compared to A,„in given in TABLE I (that is the dimensionless viscosity 
coefficient fi sufficiently small compared to unity), the effect of viscosity is also negligible. In that case, the waves are 
well described by an idealized model of a plasma being strictly cold and collisionless (Tg — — rj = 0). Within this 
model, all plasma parameters oscillate in time with the plasma frequency ujpe (1), independently of the initial density 
shape or its amplitude. Furthermore, the electron number density Tie is always greater than half its equilibrium value 
I no- Any minimum in the initial electron density distribution ne{x,t — 0) tending to |no would produce an exploding 
maximum after half a period of the plasma oscillation. These facts were known for a purely cosinusoidal deviation of 
the initial electron number density from its equilibrium value no, but are demonstrated here for arbitrary deviation. 

When increasing the wavelength A (i.e., for ly small but finite, e.g., = 0.015), growth of the maximum is reduced, 
and the minimum can fall below ^no, see Figs. 3 and 10. And for i/ relatively large, e.g., ly = 0.1, the maximum 
bifurcates as shown in Figs. 4 and 11. 

Our requirement that Ue/riQ should not be considerably larger than unity concerns the accuracy of the viscous force 
(for validity of linearisation and quasi- neutrality) . However, this term only corrects the motion, while the driving 
electric force and the convective non-linearities are taken into account exactly. In fact, only for relatively large v, 
must one stipulate that Ue should not exceed hq too much. 
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